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For translationally shape invariant potentials, the exact quantization rule proposed by Ma and 
Xu is a direct consequence of exactness of the modified WKB quantization condition proved by 
Barclay. We propose here a very direct alternative way to calculate the appropriate correction for 
the whole class of translationally shape invariant potentials. 

PACS numbers: 



In 2005, Ma and Xu |l|,|2|, on the basis of a previous work of Cao |3j,|4( and al., proposed a new improved quantization 
rule permitting to retrieve the exact spectrum of some exactly solvable quantum systems. Since, numerous papers 
have been published on the possible applications of the Ma-Xu formula to different systems }5Ml4l|. Until now, the 
link between the exactness of the Ma-Xu formula and the solvability of the considered system seems to be stayed 
unclear. In fact, the translationally shape invariance of the potential suffies to ensure the exactness of the Ma-Xu 
formula. This is a direct consequence of a result obtained first by Barclay [lj| in 1993 and largely overlooked. He 
established that for these specific potentials the higher order terms in the WKB series can be resummed, yielding to 
an energy independent correction and he showed that this Maslov index can be obtained in a closed analytical form. 
More than ten years later, Bhaduri and al. [2(]| proposed another interesting derivation of this result which rests on 
periodic orbit theory (POT) 

In the present article, after having established the connection between the Ma-Xu formula and Barclay's result, we 
propose an alternative way to calculate the Maslov index for every translationally shape invariant potential (TSIP). 
Starting from a classification of TSIP which uses new criterions [22| , equivalent to the Barclay-Maxwell ones [18(] , we 
show how to obtain this index using simple complex analysis tools. 



INTRODUCTION 



MA-XU QUANTIZATION FORMULA 



Consider the stationary one dimensional Schrodinger equation (h = 1, m = 1/2): 




where p(x) = *JE — V(x) is the classical momentum function for an energy E. If x\ and xi are the classical turning 
points then p 2 (x) > for x G [xi,^]- 
Defining: 



w{x) 



i//(x) 



(2) 



i/j(x) ' 



we have: 



w'{x) 



ijj"{x) 



+ w 2 {x). 



xjj{x) 



In other words, if ip(x) satisfies Eq.JT]), w(x) is a solution of the following Riccati equation: 



w'(x) = p 2 (x) + w 2 (x) = E - V(x) + w 2 (x). 



(3) 
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Every node of ip{ x ) (necessarily a simple zero) corresponds to a simple pole of 10(35) which decreases in the interval 
X2\- Consequently, when x runs through [11,352]? at each node of ip( x ) the function w(x) is subject to a discontinuity 
from —00 to +00. If we define the phase 9{x) via: 

tanfi(i) = -^, (4) 
w[x) 

we can write: 

( p( x ) \ 

6(x) = arctan + nir, (5) 

V W ( X )J 

where arctan y G [— 7r/2,7r/2] is the principal determination of the reciprocal of the tangent function and where n 
increases of 1 at each node of w(x). 

Then, if x\ and 22 are not nodes of tp(x): 

9'{x)dx = Ntt, (6) 



N being the total number of nodes of w(x) on [x±, X2\- 
But using Eq.(J3]) we also have: 

nil \ / \ p'(x)w(x) . 

0(x)= P (x)-^J^l. (7) 
Considering the 11 th bound state at energy E n for which N n = n + 1, we obtain the following identity: 



■r 2. 



p n (x)dx — mr + "f{E n ), (8) 



where the correction 7(£ , n ) to the lowest order WKB condition is given by the following integral: 

«E n )=n+ r n Wn{X }^ x) dx. (9) 



w' n {x) 

Ma and Xu [l], Q have observed that que for a large class of exactly solvable potentials, this integral correction is 
in fact independent of the considered energy level, and can be calculated from the ground state: 

l{ E n )= 1 (E )=, + r^)f^ d x. (10) 

Jxi,o w 0\ x ) 

Inserting Eq. tfTOl) in Eq.® gives the Ma-Xu formula: 



X2,. 



p n (x)dx = nil + j(Eq). (11) 



The knowledge of the ground state characteristics permits then to calculate the 7(^0) correction and reaches to 
an implicit formula for the other energy levels E n . The explicit calculations have been performed for many exactly 
solvable examples 5 14]. Nevertheless the question of a direct link between the exact solvability of the potential and 
the validity of the above formula is still stayed open. 

As much as we know, all the closed form analytically solvable quantum mechanical systems belongs to the set of 



translationally shape invariant potentials (TSIP) in the sense of SUSY quantum mechanics 15]. In fact for this set, the 



validity of the Ma-Xu formula follows from a result established in 1993 by Barclay [19(. In this article, he established 
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that for these potentials, the WKB series can be resummed beyond the lowest order giving an energy independent 
correction which can be absorbed in the Maslov index and written in a closed analytical form. He showed equally 
that this result is directly correlated to the exactness of the lowest order SWKB quantization condition [TBI . Hr| . The 
starting point is the definition of two classes of potentials, each characterized by a specific change of variable which 
brings the potential into a quadratic form. This two classes are shown to coincide to the Barclay-Maxwell classes [l8| 
which are based on a functional characterization of superpotentials and which cover the whole set of TSIP. In both 
cases, the action variable for an energy E n is of the form: 



In = f p n (x)dx = 2 p n (x)dx = 2 (nir + 7) , (12) 

where 7 is an energy-independent correction characteristic of the considered potential. As shown later by Bhaduri 
and al. [2(J, this results can be retrieved in an elegant way using POT [2l| . 

The translational shape invariance of the potential is then a sufficient condition for the validity of the Ma-Xu 
prescription which reaches to Eq. (|ll[) . 

In the following, we propose an alternative way to recover Barclay's result. It lies on a different characterization 
of the Barclay-Maxwell classes that we recently presented (22j. In both cases the action variable can be rewritten 
as a complex integral on an uniformization domain including only one branch-cut. Using standard tools of complex 
analysis, it can be readily calculated to recover the explicit form of 7 for the whole set of TSIP. 

Note that analogous complex analysis techni ques have been employed by Bhalla and al [l^-EBl in the frame of 
the quantum Hamilton- Jacobi (QHJ) formalism [28], [2^] where, starting from the exact QHJ quantization condition, 
they have computed the spectrum of numerous potentials for which they have also showed the exactness of the 
lowest order SWKB quantization condition. Moreover, Cherqui et al j2?J have recently shown that for algebraic and 
hyperbolic translationally shape invariant potentials, the shape invariance condition provides sufficient information 
on the singularity structure of the quantum momentum function to determine directly the energy spectrum of the 
system from the exact QHJ quantization condition. 



FIRST CATEGORY POTENTIALS 

We say that a one dimensional potential is of first category [22J if there exists a change of variable x —> u transforming 
the potential into an harmonic one V(x) — > V(u) = X2U 2 + Xiu + Ao, such that u(x) satisfy a constant coefficient 
Riccati equation: 

M x ) _ a , a ..r^ , a ..2 



dx 



A + A lU (x) + A 2 u 2 (x), (13) 



du/dx being of constant sign in all the range of values of x and u. 

The one dimensional harmonic oscillator correspond to the special case A\ = A2 = and the Morse potential is 
generically associated to the case A\ ^ 0, A2 = 0. 



Harmonic oscillator 



It is perfectly well known that the EBK quantization condition is exact for the harmonic oscillator implying then the 
exactness of the Ma-Xu formula with a Maslov index equal to 1/2. Nevertheless, we will examine this case completely 
as a first example of use of the complex variable integration technique. 

The harmonic oscillator potential with zero ground state energy is: 

V(x,u) = (cu/2) 2 x 2 - u/2. (14) 
The classical turning points ±xo, n at energy E n = ncu satisfy: 



(15) 
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The half action variable for a classical periodic orbit of energy E n is then: 



p n (x)dx^ L ^ dxJxl n -x 2 . 



(16) 



The complex function Jh{z) = \jx\ n — z 2 , defined on C— ]— Xo in ,Xo, n [t admits only one isolated singularity at 
infinity. Then: 



dxfjj(x) = 7riRes oo) . 



(17) 



The asymptotic behaviour at infinity on the Riemann sheet for which Uq n — z 2 — Xq n + y 2 > on the positive 
half imaginary axe (z = iy, y > 0) is given by: 



z 2 ~ —iz 1 

z— foo \ 



b 0,n 



■ . ix 0,n 1 1 s 

-iz+ —-- +0( — ). 
2 z z l 



(18) 



Consequently: 



Res(/ ff (z),oc) = -- 



(19) 



and: 



which is the expected result. 



p n (x)dx 



n + - I 77. 



Morse potential 



(20) 



As for the harmonic oscillator, the EBK formula is exact for the Morse pot ential. This second example is nevertheless 
very instructive. The Morse potential with zero energy ground state is [151 ]: 



V{x) =A A + B 2 e' 2ax -2B[A+ -J E - ax , a > 
Using the change of variable y — exp (—ax) , dy = —aydx, it becomes: 



(21) 



V(y) = B 2 (y- y ) 2 + V 



(22) 



with I/O = (A + a/2) /B and: 



v = a 2 -bW = (a + ^) 2 -a\ 

This is a translationally shape invariant potential with a spectrum given by [22 



(23) 



(24) 
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with a = A and a,k = A — ka. 

As for the classical turning points, they are given by: 



V% n =yo+ "o,n 
Vl,n = 3/0 - U0,n 



(25) 



with: 



B 2 

The half action variable for a classical periodic orbit of energy E n is then: 



u ,a = \ly 2 a -^- (26) 



dy 



/ p n {x)dx = - y/E n - V(y, a) 

J xi, n a Jyi,„ V 

= — /m (u) du, 

a J—Un.n 



(27) 



where: 



\l u l.n - u2 

fM(u)=V_ r . (28) 



The complex extended function Jm (z) — y u^ n — z 2 / (z + y ) T, dehncd on C— ]— u , n , uo,n[, admits two isolated 
singularities at infinity and in — y . Then: 



duf M (u) = iri (Res (f M (z) , -y ) + Res (f M (z) , oo)) . 

""0,71 

With the same choice as before for the square root determination we have: 



(29) 



z^>oo \ Z Z 

which implies: 



!m{z) ~ » (-1 + ^ + 0(1)) tm 



Res (/ M (z) , oo) = _i a + ®/ 2 _ (3i) 



The residue in — y (which is a simple pole) is readily obtained as: 



Then: 



which is the expected result. 



Res(f M (z),-y a ) = i^-^. (32) 



J p„{x)dx = + 7T. (33) 
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Effective radial potential for the Kepler-Coulomb problem 

Consider the effective radial potential for the Kepler-Coulomb. It can be written as [22j 



V(y) = 1(1 + l)y 2 - 1V + 4(/ ^ 1)2 , k > 0, (34) 



where y — l/x and dy = —ydx. We can rather use the equivalent form: 



with yo — j/2l(l + 1) and 



V(y) = 1(1 + l)(y- yof + V , (35) 



7 2 7 2 

Vo = W+W W+V' (36) 



This is a translationally shape invariant potential characterized by the following energy spectrum 221 ] : 



where a = I + 1 and a*, = £ + 1 + k,. 

The classical turning points are given by: 



V2,n = VO + U ,n 
Vl,n = VO - U0,n 



(38) 



with: 



1 / 7 2 7 2 

U ^ = 7WTT)\ WTT) 4af (39) 



The half action variable for a classical periodic orbit of energy E n is then: 



^ Pn (x)dx = - y/E n - V(y,a)% (40) 



Xl 



= -V^T+l) / ' d«/jf («) 

where: 



/«:(«) = ^ ^- (41) 

(u + yo) 



Again, the complex function fx( z ) = \J u "o n ~ z2 1 i z ~^ Vo) j defined on C— ]— Wo ;n , Uo, n [, admits two isolated 
singularities at infinity and in —yo. Then: 



ruo,n 

I duf K (u) = iri (Res (f K (z) , -y Q ) + Res (f K (z) , oo)) 

— «0,n 



(42) 
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With the same choice as before for the square root determination, we have: 



/*(*) - — (l-%)"(l + -)" 9 - ( 43 ) 

z^oo z \ Z z I \ Z / z^oo z Z z 

which implies: 



Res (f K (z) ,00) =i. (44) 
The residue in —y which is a double pole is readily obtained as: 

Bes (f K (z),-y )=i J^— ■ (45) 

Then: 

xa, „ 

p n (x)dx — mr + 7 , (46) 
where: 



which is the expected result. 



7 = 7T (l + 1 + y/l(l + 1)) 



Other first category potentials 



If we except the two preceding examples, all the other first category potentials are such that there exists a change 
of variable x — > y transforming the potential into an harmonic one: 



V(x) -)• V(y) =a{aTa)y 2 + \ t y + A (a), 



with Aq (a) = A^/4a 2 — aa (in order to have a zero energy ground state) and where y(x) satisfies a constant coefficient 



Riccati equation of the form [221 ] : 



~ =a±ay 2 (x) > 0. (47) 



The potential can still be written [221 ] : 



V±(i/,a) =a(aT«)(y-yo) 2 + 1/ , (48) 



where: 



2/o = -Ai/a (a =F a) , V = \ (a) - a (a =f a) y 2 (49) 
V is a translationally shape invariant potential and its energy spectrum is then given by [22]: 



E n (a) = (f) 1>± (a) - cj) l ± (a„) , 



(50) 



where au = a ± ka and: 



A? 



1± (a) = T a 2 + I ± 2 . (51) 
The classical turning points y^ n for an energy E n are determined by the condition: 

E n = V(y i>n ) ^ E n -fl(flTa) u 2 Qn + Vo, (52) 

where: 



V2,n = Va + UQ,n 
Vl,n =U0- U Q , n 



(53) 



and: 



^.n-WTi + yg- Vi° n v (54) 

Consider first the case where dy/dx = a + ay 2 (x). The half action variable for a classical periodic orbit of energy 
E n is given by: 



X2. n /•Ste.ra 



p„(x)efa; = / v/£; n - V(y,a) — — — - (55) 

a{y-i)(y + i) 



Xl, n J VX,n 

Uo,n 

fi,+(u)du, 



y/a(a — a) 



a 

where: 



«0, 



t 2 _ z 2 

^- (.:Jc-o (56) 

with zq — —yo + i. 

The complex extended integrand /i + (z), defined on C— ]— uo, n , uo,n[, admits three isolated singularities at infinity, 
Zo and Zo- Then: 

"0,71 

dufi,+ (u) = ni (Rcs(/i i+ (z),z ) + Res (/i,+ (z), 2o) + Res oo)) . (57) 

-U0,n 

When z — > oo: 

^ -*+0(^), (58) 

which implies: 

Res(/i,+ (z),oo)=i (59) 
As for the residues in z and Zo, using Eq.([5"2"]). Eq.(J53J), Ea.([5"l|) and Eq. lJSD]) we deduce: 

Res(/i,+ (z)^o) = -— 7 = A ian + ^~) ( 60 ) 
2 v / a(a - a) \ 2a„ / 
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and: 



Then: 



Res(/i, + O),z ) = - 1 = \-ia n + ^- ) . (<>1 ) 

{a -a) \ 2a„ 



X2, n 

p n (x)dx = nir + 7 (62) 



with: 



a 

,2 



7ra / / a , 

7 = — ■ (63) 



Consider now the case dy/dx = a — ay (x). The half action variable for a classical periodic orbit of energy E n is: 



%2, n PV2 



p n (x)dx = I y/E n - V(y,a)— ^— — - (64) 

Jyi, n a{l-y){l + y) 

yja (a - a) f u °-' 



a 

where: 



fi-(u)du, 



U 0,n - u2 



= , w V (65) 

(u - zi) (u - z 2 ) 

with z\ = — 1 — z/o an d ^2 = 1 — t/o • 

The complex extended integrand /^(z), defined on C— ]— uo,n, wo,n[> admits three isolated singularities at infinity, 
z\ and Z2. Then: 

"0,7! 

dufi-(u) = iri (Res (/1 -(z),Zi) + Res (/1 _(z), z 2 ) + Res (/1 _(z), 00)) . (66) 

-M0,n 

When z — >• 00: 

/!,+(«) ^ - 1 +0(^), (67) 

which implies: 

Res 00) =t. (68) 

As for the residues in zo and zq, from Ea. (f52"|) . Ea. (|9"Tj) . Ea. (l92|) and Eg. ((50)) we deduce: 

Res (fi,-(z), z 2 ) = % (a n + £-) (69) 
2y/a(a - a) \ 2a„J 

and: 
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Note that in this case, we had to change the square root determination when we pass from z\ to z%, since we have 
either -1 < y 1>n < y 2 , n < l(when y = tanh(ax + ^)), that is -1— y < —uo, n < Uo,n < 1 — 2/0, or j/ 1)W < -1 < 1 < y 2>n 
(when y = coth(ax + <f)), that is — ito,n < — 1 — yo < I — yo < «o,n- 

Then: 

X2, » 

p n (x)dx = mr + 7 (71) 

Xl, n 

with: 



a \ v a 

We recover the results obtained in [2(| for the first class potentials. 



SECOND CATEGORY POTENTIALS 



ira I / a . 
7 = 1 + 4/1- - . (72) 



We say that a one dimensional potential is of second^ category if there exists a change of variable x — > u transforming 
the potential into an isotonic one V(x) — > V(u) — X 2 u 2 + Ao + ^f, such that u(x) satisfies a constant coefficients 
Riccati equation of the form: 

*P.=A +A 2 u*(x), (73) 
ax 

du{x)/dx being of constant sign in all the range of values of x and u. 

Isotonic oscillator 



The isotonic potential with a zero energy ground state (E Q = 0) is [22j, [30j: 



1 "f-'-'V' (74) 



Its energy spectrum is given by [221 ] : 

£ n (a) = 2nw, (75) 

where a = (f ,Z + l) and a k = (f ,1 + 1 + fc). 

The classical turning points iTj in for an energy i?„ are determined by the condition: 

E n = V(x z . n ) «■ - w ^2n + i + x?„ + 1(1 + 1) - 0, (76) 

that is: 



x 2,n — u 0,n + &n 



(77) 



with: 



2n + 1 + - 2// 3\ 2 

uo,« = 2 5 n = - W 2n + Z + - - Z(Z + 1). (78) 



J /* l/l O 
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The half action variable for a classical periodic orbit of energy E n is then: 



p n (x)dx = — I duf^u), 



where we have defined u — x 2 and: 

J( u - x l.n) { x l,n - U ) 

fi(u) = 



The complex extended function fi(z) = y (z — x\ „) {x\ n — z) /z, dehned on C— ] x\ n , x\ n [, admits two isolated 
singularities at infinity and in 0. Then: 



dufjiu) = iri (Res (//(»), 0) + Res (//(z), oo)) . (79) 
When z — >• oo: 

fiiz) - -i (l- Xln + Xln +0(\)), (80) 

which implies: 

2n + I + | 

Res (fi(z), oo) = -iu ,„ = -2i = — ±, (81) 



As for the residue in 0, from Eg. ([5^)) we have: 



Then: 



Res (/z(*),0) = J-x 2 1>n xl n = i-VW+Tj. (82) 



p„(x)dx = mr + 7, (83) 

Xl.n 



with: 



(84) 



Other second category potentials 



If we except the preceding example, all the other second category potentials are such that there exists a change 
of variable x — »■ y transforming the potential into an isotonic one V(x) — > V(y) — A2?/ 2 + p- + \q, where y{x) > 
satisfies a constant coefficient Riccati equation of the form [22]: 

^=a±ay 2 (x), (85) 
dx 



dy{x)/dx being of constant sign in all the range of values of x and y. 
As shown in [221 ] . the potential can be written: 



where a — (A, fi) and: 



V(y, a) = A (A T a) y 2 + M(M 2 " } + A (a) , 



Aq (a) = —a (A ± fx) — 2A/i. 



V is a translationally shape invariant potential and its energy spectrum is given by [22 



E n (a) = ± {<j) 2 ± (a n ) - (j) 2 ± (a)) , 
where a = ao, at = (A&, /%) = (A ± fea, /i + fca) and: 

02,± ( a ) = 02,±(^ M) = ( A ± M) 2 ■ 

The classical turning points y^ n for an energy E n are determined by the condition: 

E n = V(vi, n ) «A(At«) yf in + (A (a) - £„) y\ n + p (fi - a) = 0, 

that is: 

f 2/2,n = u 0,« + 5 n 

with: 



_ g„ - A (a) _ U E n - \ (a) \ 2 fi(fi-a) 
UQ < n ~ 2A(At«) ' V V 2A(At«) J A (A =F a) ' 

In the case dy/dx = a + ay 2 , the half action variable for a classical periodic orbit of energy E n 



p n (x)dx = / \/E n - V(y,a) 

J 2/1,71 

(A -a) p 2 



r/i/ 



2a 



duf 2 ,+ {u), 

Vl.n 



where we have defined u — y 2 and: 



/ 2 ,+ (*) 



z(z + l) 

This last integral is readily calculated by using a complex variable formalism. We have indeed: 



»2,n 

duf 2 , + {u) = iri (Res (/ 2 ,+ («),0) + Res (/ 2 ,+(«), -1) + Res (/ 2 , + (z), oo)) 
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When z — > oo : 



/ a ,+(*) ^ (96) 



z^oo z 'z 2 ' 



which implies: 

Res(/ 2 ,+ (z),oc) =i (97) 
As for the residues in and —1, from Ea. (j9"U]) we have: 



Res (£,+(*), 0) = ^-vl n yl n = iJ x{x_ay ( 98 ) 



Using Eq.(j90]), Eq.(JMJ», Eq.(|92|) and Eq.((88j), we deduce (-1 < < y\ „ < j/f J: 



Then: 



Res (/ 2 , + (z), -1) = - v/- (l + 3/?,„) (j/ 2 2 n + l) = %=^ ' (99) 



p n {x)dx = nir + 7 (100) 

2:1,71 



with: 



A similar analysis can be led in the case dy/dx = a — ay 2 (x). If we except the specific case of the Scarf II potential 



22j . we can write: 



X2 ' n p n (x)dx = - ^ A ( A + Q) T 2 '" du/ a ,_(t*), (102) 



where we have defined u = y 2 and 



We have: 



When z — >• 00: 



/».-(*) = ^7" T) ' ( 103 ) 



y 2 ,n 

duf 2 ,-(u) = Tri (Res (/ 2 ._ (z), 0) + Res (/ 2 ,_(z), 1) + Res (/ 2 ,_ (z), oo)) . (104) 



z->oo z Z z 
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which implies: 

Res(/ 2 ,_(z),oo) =%. (106) 
As for the residues in and 1, using Eq.®, Eq.([HI]), Eq.® and Eq.([55]), we deduce: 



Res (/ 2> _(*),0) = -^-Vl n v\ n = -J fj^j (107) 



and 



Ree (&_(*), 1) = y/(l-v 2 1 , n )(ti, n -l) = %=4 - (108) 
v v /A(A + a) 

Note that in this case, we had to change the square root determination when we pass from to 1, since we have 
< yl. n < yln < 1 ( V = tanh(air + if)). 
Then: 

X2,n 

p n (x)dx — nn + 7 (109) 
with: 



In the case of the Scarf II potential, the variable y(x) = tanh (ax/2 + itt/A) is a pure phase factor (\y\ = 1) and the 
path of integration in the action variable surrounds the branch cut which is an arc of the unit circle. Nevertheless, 
we can follow the same reasoning as before and we recover the result given in Eq. ()110p . 

Eq. (|10ip and Eci. (|110[) correspond to the results obtained in [20( for the second class potentials. 

CONCLUSION 

We have shown how to calculate exactly and in a general way the action variable for the whole set of translationally 
shape invariant potentials. The correction to lowest order WKB formula appears as a constant term redefining the 
Maslov index. This result implies immediately the exactness of the Ma-Xu formula for every TSIP. The employed 
techniques of complex analysis are standard but even so instructive. The two basics cases are the harmonic and 
isotonic ones (for which we have an equispaced quantum spectrum and isochronicity at the classical level) . In this two 
cases, the action integral involved is a sum of at most two residues and the linear n dependence (n being the energy 
quantum number) comes from the residue at infinity. For all the other cases, the integral involved is obtained from 
the two basic cases by deforming the integration measure with a weight which is the inverse of an at most quadratic 
polynomial. The residue at infinity becomes energy independent and the linear n dependence takes its origin in the 
finite poles (with an eventual compensation of nonlinear n dependent contributions) . 
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